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ABSTRACT 


The genus Macrocarpaea (Griseb.) Gilg (Gentianaceae, Helieae) is among the largest woody genera of tropical gentians, with 
most of its species occurring in the wet mountainous forests of the Andes. Phylogenetic and dispersal-vicariance analyses 
(DIVA) of 57 of the 105 currently recognized species in the genus, using two data sets from nuclear DNA (ITS and 5S-NTS 
sequences) and morphology. show a single origin of the Andean species from an ancestral distribution that includes 
southeastern Brazil. Within the Andes, species divide into two major clades: (1) northern species from the cordilleras of 
northern Ecuador, Colombia, and Venezuela; and (2) southern species of the Andean Amotape-Huancabamba Zone in 
Ecuador and Peru, as well as the Andes of central and southern. Peru and Bolivia. The Amotape-Huancabamba Zone is 
supported as the ancestral area for Macrocarpaea within the Andes. There are repeated speciation patterns within the Andes, 
and three Mesoamerican species derive from the northern clade, as is the single sampled species from the Guayana Shield. 
The position of the subclade of the three Caribbean species is less certain, but it currently nests among Andean species. An 
Atlantic coastal Brazilian clade is placed as sister group to all other Macrocarpaea, providing further support for an ancestral 
refuge in southeastern Brazil for the Helieae. The biogeographic analysis showed that local speciation is more common than 
long-distance dispersal, and allopatric speciation is more common than sympatric speciation. Using detailed, georeferenced 
herbarium collection data, patterns in environmental characteristics between clades and sister species were analyzed with 
Spatial Evolutionary and Ecological Vicariance Analysis (SEEVA), utilizing geographic information system (GIS) and 
statistical methods. Sister clades and taxa were evaluated for statistical significance in variables such as annual rainfall and 
temperature, elevation, temperature and rainfall seasonality, geological bedrock age. and soil type to evaluate ecological 
vicariance between sister groups. The results indicate that there are no general patterns for each variable, but that there are 
many significant divergences in ecological niches between both larger sister groups and sister species, and ecological niche 
conservation was also observed when subsequent nodes in the phylogeny were compared. 

Key words: Biogeography, ecology. Gentianaceae, Macrocarpaea, Neotropics, niche, South America, speciation, 
vicariance. 


The tropical Andes are one of the most biologically due to recent geological uplift. Biodiversity in South 
diverse areas in the world and a biodiversity hotspot America and Andean biogeography has been reviewed 
(Myers et al., 2000; Rodriguez-Mahecha et al., 2004a, by Burnham and Graham (1999) and Young et al. 
b; Jorgensen et al., 2007; Morawetz & Raedig, 2007), (2002), highlighting the complex history, topography, 
despite being a relatively young part of South America and vegetation. patterns of the Andes. The high 
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biodiversity has been explained through habitat 
diversity resulting from large differences in bedrock 
types, soils, climate, and elevation, which are also 
fragmented due to the dissected topography. Species 
have spread north and south along this jagged 
mountain chain, and also up and down in elevation 
during colder times through repeated glaciations. Both 
allopatric and sympatric speciation processes could 
have resulted in such large species diversity. 

Generally speaking, the theory behind sympatric 

speciation suggests that competition among popula- 

tions led to subsequent ecological niche divergence, 
whereas allopatric speciation invoked a spatial barrier 
and not necessarily any ecological niche separation 

(Mayr, 1963). This is applicable not only to extant 

sister species pairs, but also to deeper clades in the 

phylogeny that represent ancient allopatric or sym- 
patric events. 

Progress in methodological developments has now 
made it possible to evaluate phylogenetic patterns and 
speciation processes in a historical and spatial context 
including ecological data. In this study, we use 
molecular and morphological data for phylogenetic 
reconstruction using parsimony criteria, dispersal- 
vicariance analysis (DIVA) for reconstruction of 
ancestral areas (Ronquist, 1997), and Spatial Evolu- 
tionary and Ecological Vicariance Analysis (SEEVA) 
(Struwe, 2008). 

In this paper, we will evaluate several hypotheses 
based on a comprehensive data set from Macrocarpaea 
(Griseb.) Gilg. 

l. Continental patterns: What is the relationship 
between species in the Andes to species in 
Mesoamerica, the Caribbean, southeastern Bra- 
zil, and tepuis of the Guayana Shield? 

2. Andean patterns: What is the ancestral distri- 
bution area in the Andes and what large-scale 
Andean biogeographic patterns are present? 

3. Speciation patterns: Is sympatric or allopatric 
speciation most common in Andean species? Is 
dispersal or vicariance most common in Andean 
species? 

4. Ecological niche patterns: Does allopatric spe- 
ciation indicate a relatively larger shift in 
ecological niche characteristics than sympatric 
speciation that occurs within the same area? Do 
vicariant (allopatric) versus overlapping (sym- 
patric) sister taxa have similar or divergent 
ecological niches? 

Our study group, Macrocarpaea (Gentianaceae, 
Helieae), occurs in mountainous regions of the 
Neotropics at (30—)1500-3000(-3800) m elevation. 
It is composed of woody shrubs, small trees (up to 
10 m), and herbs (one species; M. rubra Malme) with 
large (up to 7 em long), campanulate, white, yellow, to 


green flowers that are visited diurnally by humming- 
birds, butterflies, and insects, and nocturnally by bats 
and moths (Grant & Struwe, 2001; Grant, pers. obs.). 
Its woody habit is uncommon in the Gentianaceae and 
led to a comparative study of its wood anatomy 
indicating secondary derivation of woodiness from 
herbaceous ancestors (Carlquist & Grant, 2005). With 
105 species currently recognized, Macrocarpaea is by 
far the most species-rich genus in the Helieae, while 
the majority of genera in this tribe are either 
monotypic or have less than 10 species, e.g., 
Celiantha Maguire and Yanomamua J. R. Grant, P. 
J. M. Maas & Struwe (Struwe et al., 2002; Grant et al., 
2006). During monographie studies on Macrocarpaea, 
more than 3500 herbarium sheets were examined, 
resulting in the description of 70 new species (Grant 
& Struwe, 2001, 2003; Grant, 2003, 2004, 2005, 
2007, 2008; Grant & Weaver, 2003). Fieldwork by 
Grant, especially in the Huancabamba region of 
Ecuador and Peru, continues to uncover new species 
that are often cryptically different from one another 
when seen as herbarium specimens, but are otherwise 
clearly distinct in the field. Likewise, molecular 
phylogenetics of Macrocarpaea (Struwe et al., 2009) 
has given strong support for its monophyly, for the 
establishment of an infrageneric classification (Grant, 
2005), and to the understanding of difficult species 
complexes. 

Macrocarpaea has a broad distribution in moun- 
tainous regions of the Neotropics and is comprised of 
generally narrowly endemic species. It occurs princi- 
pally in the Andes (87 species from Venezuela, 
Colombia, Ecuador, Peru, and Bolivia), with outlying 
species in the Pantepui of the Guayana Shield (six 
species), in southern Mesoamerica (five species in 
Costa Rica and Panama), in the Atlantic forest of 
southeastern Brazil (four species comprising section 
Tabacifoliae Ewan), and in the Greater Antilles of the 
Caribbean (three species: one each from Cuba, 
Hispaniola, and Jamaica). The three main centers of 
diversity are in the Colombian Andes (section 
Macrocarpaea, where 58% of the species of that 
section occur), the Amotape—Huancabamba region 
straddling southern Ecuador and northern Peru 
(section Choriophylla (Griseb.) J. R. Grant, 74% of 
its species), and the central Peruvian Andes (section 
Magnolifoliae Ewan, 50% of its species). The genus is 
absent from lowland Amazonia, Mexico, Mesoamerica 
north of Costa Rica, and southern temperate South 
America. Here, we use Macrocarpaea as a model for 
understanding South American, and particularly 
Andean, speciation, biogeography, and dispersal. 
Because Macrocarpaea has many narrowly distributed 
Andean endemics, it is an excellent candidate for this 


kind of study. 
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MATERIALS AND METHODS 
PHYLOGENETIC DATA AND ANALYSES 


Fifty-seven of 105 total species (54%) of Macro- 
carpaea, representing all sections as defined by Grant 
(2005), were included in the phylogenetic analysis 
(Table 1). Five outgroups were selected for orientation 
of the phylogenetic trees (Nixon & Carpenter, 1993), 
based on previous phylogenetic work in the Gentia- 
naceae and all of tribe Helieae using trnL and matK 
chloroplast data as well as nuclear ITS markers 
(Struwe et al., 2002, 2009). Outgroups used include 
Symbolanthus jasonii J. E. Molina & Struwe from the 
Symbolanthus G. Don subclade of Helieae, and 
Chorisepalum psychotrioides Ewan and three species 
of Tachia Aubl. (T. grandiflora Maguire & Weaver, T. 
guianensis Aubl., and T. occidentalis Maguire $ 
Weaver) assigned to the Macrocarpaea subclade 
(Struwe et al., 2002). 

The DNA data in this study were developed for a 
detailed phylogenetic study of the genus by Grant. For 
molecular data we used two different nuclear genomic 
areas, ITS and NTS (for 55 RNA [5S-NTS]). These 
regions successfully resolved phylogenetic 
relationships in tribe Helieae on both generic and 
tribal scales in the past (Struwe et al., 2002, 2009; 
Gould & Struwe, 2004). Additional sequences for 57 
species of Macrocarpaea and the selected outgroups 
were downloaded from GenBank (see Table 1 for 
GenBank numbers). In total, 47 ITS and 56 5S-NTS 


taxon-specific sequences were obtained from Macro- 


have 


carpaea. Sequences were aligned using ClustalX 
(Thompson et al., 1997) and Se-Al v.2.0 for Macintosh 
(Rambaut, 2002). 

In addition to molecular data, nine morphological 
characters were coded for all included species 
(Table 2), and these were primarily selected from 
pollen and seed surface structure (Grant & Struwe, 
2001, 2003; Grant, 2003, 2004, 2005, 2007; Grant & 
Weaver, 2003). Four morphological seed types were 
identified by Grant (2005), including two types 
without wings (flattened and rimmed), and two with 
wings (perimetrically winged and winged). Character 
l refers to the general shape of the seeds, varying from 
three-dimensional (rectangular to spheroid) to com- 
pletely flattened, as does character 4 separating seed 
types with either a square or long-linear outline. 
Character 2 refers to the presence or absence of seed 
wings, and character 3 is coded based on whether 
these wings occur on opposite ends as in a bowtie 
(winged type) or surround the seed entirely (perimet- 
rically winged). Character 5 is the weight of an 
individual seed in milligrams, generally less than 
0.1 mg in Macrocarpaea and greater than 0.1 mg in 


Chorisepalum Gleason & Wodehouse, Tachia, and 


Symbolanthus. Character 6 refers to features of the 
pollen exine surface, and character 9 identifies 
whether the pollen is shed as monads or as tetrads 
or polyads. Character 7 differentiates the axillary 
flower position of Tachia, as compared to all other 
genera that have multiflowered terminal inflorescenc- 
es, and character 8 codes for the color of the corollas. 
Character states were treated as unordered, and 
characters were regarded of equal weight. The two 
DNA alignments were combined with the data matrix 
of nine morphological characters and analyzed 
together simultaneously. 

The phylogenetic, unconstrained parsimony analy- 
sis was performed in Winclada and NONA (Goloboff, 
1999; Nixon, 1999, 2002) using 500 random repli- 
cates, five starting trees per replicate, MaxTrees as 
10,000, and the tree-swapping algorithm was multiple 
tree bisection-reconnection (TBR) + TBR. Bootstrap- 
ping was done with 300 replicates in PAUP* v.4 
(Swofford, 2000). 


COLLECTION DATA 


A total of 794 georeferenced herbarium collection 
records of the 57 Macrocarpaea species included were 
entered into an existing FileMaker Pro 7.0 (Claris Pro, 
Santa Clara, California) database. Each locality was 
georeferenced using printed and online maps, atlases, 
and gazetteers. Identification of all collections was 
confirmed by Grant, and an exsiccatae list with the 
georeferenced coordinates is available upon request. 
Only recorded locations known to the nearest minute, 
nearest second, or label with global positioning system 
(GPS) coordinates were included in the analysis. 
Coordinates were converted to a point shapefile in 
ArcGIS v. 9.2 (ESRI, Redlands, California), and 
distribution maps of each species and the genus were 
produced. 


BIOGEOGRAPHIC ANALYSIS 


Nodes present in the selected most parsimonious 
tree were classified as either representing sympatric, 
partially sympatric, or allopatric species distributions 
based on data from individual species maps. This 
classification was used in comparing impact indices 
and P values for environmental species and clade data 
during the SEEVA analysis. 

For the DIVA analysis, the total genus distribution 
was divided into 10 biogeographic units based on age, 
spatial structure and geographic connectivity, and 
geological coherence. The identified areas were: A, 
southeastern Brazil; B, Pantepui of the Guayana 
Shield; C, Greater Antilles of the Caribbean; D, 


Mesoamerica; E, Cordillera Central in Colombia and 
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Cordillera Oriental in Ecuador; F, Chocó and 
Cordillera Occidental in Colombia and Ecuador; C, 
Cordillera Oriental and Mérida in Colombia and 
Venezuela; H. Amotape-Huancabamba region in 
Ecuador and Peru; I, Cordillera Central in central 
Peru; J, Bolivia and Cordillera Central in southern 
Peru; K, Amazon Basin in Brazil, Colombia, Peru, and 
Venezuela. À map of the areas in northwestern South 
America is presented in Figure 1. 

Areas were delimited based on spatial, historical, 
and geological features; some outliers were clearly 
disjunct, such as the Caribbean (C), tepuis on the 
Guayana Shield (B), and southeastern Brazil (A). For 
areas that are part of the Andean mountain chain and 
the Mesoamerican land bridge, their distinction is less 
obvious. The boundary between Mesoamerica (D) and 
the Chocó and Cordillera Occidental in northwestern 
Colombia (F) follows evidence that lowland Darién in 
eastern. Panama has a stronger geological and 
historical connection to northern South America than 
to western Panama (Clapperton, 1993). The three 
cordilleras in Colombia and two in northern Ecuador 
were treated as separate areas because deep valleys 
separate the mountain ranges. The Cordillera de 
Mérida range in Venezuela is an extension of the 
Cordillera Oriental of Colombia and is included in the 
same area. The division between the Andean area of 
the Amotape-Huancabamba Zone (H) to the south and 
the Andean areas in northern Ecuador and Colombia 
(E + F) is correlated with the Amotape Cross, a 
geological shear zone at ca. 3^8 in the continental 
crust of the South American Plate (Clapperton, 1993). 
The southernmost limit of area H is from Río Chicama 
on the western side and Río Huallaga on the eastern 
side of the Andes and corresponds to a disjunction 
zone found in other biogeographic studies (Weigend, 
2002). The two southernmost Andean areas (I, J) in 
Peru and Bolivia are separated at the latitude of Río 
Pisco and Río Entero, which corresponds to the 
geological Abancay Deflection, another shear zone 
(Clapperton, 1993). 

A dispersal-vicariance analysis was run using the 
software DIVA v. 1.1 (Ronquist, 1996), limiting the 
maximum number of optimized areas at each node to 
three (maxareas — 3) and optimizing dispersal and 
vicariance events on a selected most parsimonious 
tree, because DIVA requires a fully resolved tree for 
its analysis. This tree was selected from all most 
parsimonious trees by being overall the most similar to 
the majority-rule consensus tree. The DIVA distribu- 
tion matrix is shown in Table 3. The DIVA optimi- 
zation method uses the parsimony criterion to limit the 
numbers of events (— steps) per tree based on the 
distribution of each species. Ancestral areas were 
optimized at each internal node in a way that limits 


vicariance and dispersal events as far as possible (i.e., 
it provides the simplest explanation of the current 
distribution data given the phylogeny). At several 
nodes, there were several equally likely area optimi- 
zalions and we selected from among these based on 
additional data from geological and spatial informa- 
tion. 


SEEVA ANALYSIS 


The theoretical background and statistical expla- 
nations and justifications for SEEVA are outlined in 
Struwe (2008) and Heiberg and Struwe (2008), with a 
short overview given here. Environmental data were 
extracted using GIS from all georeferenced collection 
localities of each species. Each variable was divided 
into four or five categories that were either quantita- 
live (e.g., rainfall amounts, in equally sized quartiles 
for the data overall) or qualitative (e.g., soil types) 
depending on data set. A table was prepared in 
Microsoft Excel (Microsoft, Redmond, Washington) 
that lists the number of collections per species for 
each variable (e.g., annual rainfall) and within each 
variable, for each category (e.g., 1-340, 341-732, 
733-1277, 1278-4000 mm). The table was imported 
into the software SEEVA v. 0.4 (Heiberg, 2008) for 
statistical analysis. The selected phylogenetic tree 
was then imported into SEEVA. 

Environmental point data for nine variables were 
extracted for the 794 Macrocarpaea collections using 
the following base layers in ArcGIS 9.2: elevation 
(U.S. Geological Survey; unit: m; format: grid; scale: 
30 are second), soil type (ESRI/ArcAtlas; format: 
vector; scale: 1:5-10,000,000), geology (unit: bedrock 
geologic age; U.S. Geological Survey; format: vector; 
scale: 1:1—5,000,000), annual mean temperature 
(BIOL; unit: ^C X 10; format: grid; scale: 1 km), 
temperature seasonality (BIO4; standard deviation X 
100; format: grid; scale: 1 km), minimum temperature 
of coldest month (BIO6; unit: ^C X 10; format: grid; 
scale: 1 km), annual precipitation (BIO12; unit: mm; 
format: grid; scale: 1 km), precipitation of driest 
month (BIO14; unit: mm; format: grid; scale: 1 km), 
precipitation seasonality (BIO15; unit: mm; format: 
grid; scale: | km), and precipitation seasonality 
(coefficient of variation; format: grid; scale: 1 km). 
All climate data were retrieved from WORLDCLIM 
(Hijmans et al., 2005). 

For each node in the phylogenetic tree, extracted 
environmental data were pooled to represent data for 
each monophyletic group (clade) and compared with 
their sister group for each node. Ecologic data are 
measured for statistically significant differences be- 
tween clades based on the distribution in the four or 
five categories. All variables were analyzed indepen- 
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Table 1. Voucher and GenBank accession numbers for 55-NTS and ITS sequences of Macrocarpaea and outgroups used 


for the phylogenetic analysis. N/A indicates sequences not available. 


Taxon Voucher 5S-NTS ITS 
Macrocarpaea angelliae J. R. Grant € Ecuador, J. R. Grant 4289 (NY) EU541681 AY397760, 
Struwe AY397761 

M. apparata J. R. Grant & Struwe Ecuador, J. R. Grant 4002 (NY) EU541683 — DQ401413 
M. arborescens Gilg Ecuador, J. R. Grant 4084 (NY) EU541686 | EU528076 
M. auriculata Weaver & J. R. Grant Costa Rica, R. L. Wilbur & Almeda 16828 (F) EU541688 N/A 

M. bangiana Gilg Bolivia, S. G. Beck 8745 (NEU) EU541690 — EU528078 
M. bubops J. R. Grant & Struwe Ecuador, J. R. Grant 4046 (NY) EU54.1692 EU528081 
M. chthonotropa J. R. Grant Peru, V. Quipuscoa 2044 (F) EU541694 N/A 

M. cinchonifolia (Gilg) Weaver Bolivia, S. G. Beck 24780 (NY) EU541696 — EU528084 
M. cochabambensis Gilg-Ben. Bolivia, A. Gentry 44200 (NY) EU541697  EU528085 
M. densiflora (Benth.) Ewan Colombia, K. von Sneidern 2523 (S) EU541700 — EU528087 
M. dies-viridis J. R. Grant Ecuador, J. R. Grant 4352 (NY) EU541702 — EU528089 
M. domingensis Urb. & Ekman Dominican Republic. D. Kolterman s.n. (JBSD)  EU541704 — EU528091 
M. elix J. R. Grant Ecuador, G. Harling & Andersson 23442 (MO) EU541706 | EU528094 
M. ericii J. R. Grant Peru, E. Rodriquez 2926 (HUT) EU541707  EU528093 
M. fortisiana J. R. Grant Peru, D. McCarroll 128 (NY) EU541709 — EU528095 
M. gattaca J. R. Grant Ecuador. J. R. Grant 4209 (NY) EU541710 — DQ401414 
M. gaudialis J. R. Grant Colombia, R. E. Weaver 2650 (GH) EU541713 — EU528097 
M. glabra (L. f.) Gilg Colombia, J. R. Grant 4310 (NY) EU541714  EU528098 
M. glaziovii Gilg Brazil. B. Rezende Silva 1318 (NEU) EU541774 = N/A 

M. gondoloides J. R. Grant Ecuador, G. Tipaz 1051 (MO) EU541716 N/A 

M. harlingii J. S. Pringle Ecuador. J. R. Grant 4049 (NY) EU541721 — EU528104 
M. innarrabilis J. R. Grant Ecuador. F. Luisier 2 (LOJA) EU541723  EU528106 
M. jactans J. R. Grant Ecuador. J. Clark 8927 (NY) EU541725  EU528108 
M. jensii J. R. Grant & Struwe Ecuador, J. R. Grant 4047 (NY) EU541724  EU528107 
M. kuelap J. R. Grant Peru, J. R. Grant 3942 (NY) EU541726 — EU528109 
M. lenae J. R. Grant Ecuador. J. R. Grant 4013 (NY) EU541727  EU528110 
M. luna-gentiana J. R. Grant & Struwe Ecuador, J. R. Grant 4028 (NY) EU541728  EU528111 
M. luteynii J. R. Grant & Struwe Colombia, 1. Cabrera & H. van der Werff EU541730 — N/A 

15769 (U) 

M. macrophylla (Kunth) Cilg Colombia, J. R. Grant 4312 (NY) EU541735 — EU528113 
M. maguirei Weaver & J. R. Grant Peru, B. Maguire 61569 (NY) EU541736 — EU528114 
M. micrantha Gilg Peru, J. R. Grant 3966 (NY) EU541737  EU528116 
M. neblinae Maguire & Steyerm. Venezuela, B. Maguire 36886 (NY) EU541739 N/A 

M. nicotianifolia Weaver & J. R. Grant Colombia, A. S. Barclay 3402 (US) EU541740 — EU528118 
M. noctiluca J. R. Grant & Struwe Ecuador, J. R. Grant 4003 (NY) EU541742 EU528121 
M. normae J. R. Grant Peru, K. García 267 (NEU) EU541743 — EU528122 
M. obtusifolia (Griseb.) Gilg Brazil. W. Thomas 14304 (NY) EU541775  EU528125 
M. opulenta J. R. Grant Ecuador, J. R. Grant 4347 (NY) EU541746 EU528124 
M. ostentans J. R. Grant Peru, B. Wallnófer 12968 (U) EU541747  EU528128 
M. pachyphylla Gilg Colombia, M. L. Bristol 1429 (GH) EU541748 — EU528129 
M. pachystyla Gilg Peru, J. Schunke-Vigo 5298 (NY) EU541751 — EU528130 
M. pajonalis J. R. Grant Peru, M. Weigend 545 (NY) EU541749 N/A 

M. papillosa Weaver & J. R. Grant Venezuela, R. E. Weaver 2629 (GH) EU541750 — EU528131 
M. pinetorum Alain Cuba. Bisse 49708 (HAJB) EU541753 EU528133 
M. pringleana J. R. Grant Ecuador. F. Luisier 1 (NY) EU541755  EU528134 
M. revoluta (Ruiz & Pav.) Gilg Peru, M. Weigend 5288 (NY) EU5417560 | EU528135 
M. robin-fosteri J. R. Grant Peru, M. Weigend 5777 (B) EU541757  EU528136 
M. rubra Malme Brazil. M. Reginato 755 (NEU) N/A EU528138 
M. sodiroana Gilg Ecuador, J. R. Grant 4210 (NY) EU541758 — EU528140 
M. stenophylla Gilg Peru, J. R. Grant 3932 (NY) EU541762  EU528142 
M. subcaudata Ewan Costa Rica, R. L. Wilbur & Almeda 16828 EU541763 = EU528143 

(DUKE) 

M. subsessilis Weaver & J. R. Grant Ecuador, J. R. Grant 4020 (NY) EU541765 EU528144 
M. tahuantinsuyuana J. R. Grant Peru, F. Woytkowski 6672A (MO) EU541766 N/A 
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Table 1. Continued. 

Taxon Voucher 5S-NTS ITS 

M. thamnoides (Griseb.) Gilg Jamaica, P. Acevedo-Rodríguez 9700 (NY) EU541767 N/A 
M. valerioi Standl. Costa Rica, S. Hill 17751 (NY) EU541769 X EU528148 
M. viscosa (Ruiz & Pav.) Cilg Peru, M. Weigend 5429 (NY) EU541770 — EU528149 
M. weigendiorum J. R. Grant Peru, M. Weigend 5363 (NY) EU541771 EU528150 
M. zophoflora Weaver & J. R. Grant Peru, J. J. Wurdack 1618 (NY) EU541772 N/A 
Chorisepalum psychotrioides Ewan N/A EU709793 
Svmbolanthus jasonii J. E. Molina Ecuador, J. R. Grant 4350 (NY) N/A EU528151 

& Struwe 
Tachia grandifolia Maguire & Weaver DQ401429 — DQ401418 
T. guianensis Aubl. DQ401430 — DQ401420 
T. occidentalis Maguire & Weaver DQ401427 | DQ401423 


dently through a modified chi-square analysis that 
gives an impact index (i) that measures the magnitude 
in trends and differences between clades on a scale 
from O to 1 for each variable and each node. P values 
were calculated using Fisher's exact test. Statistical 
significance was established at P « 0.05. 


Rest LTS 
PHY LOGENETIC DATA AND ANALYSES 


The lengths of the alignments were 687 nucleotides 
for ITS and 599 nucleotides for 5S-NTS. The 
combined molecular analysis yielded 160 most 
parsimonious trees (1252 steps, consistency index 
[CI] — 0.65, retention index [RI] — 0.74; Fig. 2). Five 
major clades are identified as monophyletic groups 
(Fig. 2, I-VI). Clade I represents three species with 
flattened seeds and pollen of the Glabra-type from 
southeastern Brazil (Macrocarpaea sect. Tabacifoliae 
sensu Grant, 2005; for distribution data, see Figs. 3, 
4) and is sister to the rest of the genus. Next is a large 
dichotomy that is poorly supported between northern 
species with rimmed or perimetrically winged seeds 
and Glabra-type pollen (clades Il and III) and 
southern species with winged or perimetrically winged 
seeds and Corymbosa- or Glabra-type pollen (clades V 
and VI). The Caribbean clade (clade IV) with rimmed 
seeds and Glabra-type pollen is weakly supported as 
sister to the clade with the more southern species. 
Clades II and HI are each monophyletic, potentially 
sister groups, and include species primarily from the 
northern Andes, the tepuis, and Mesoamerica. Clades 
II, III, and IV correspond together to Macrocarpaea 
sect. Macrocarpaea sensu Grant (2005). The only 
sampled species of six from the tepuis, M. neblinae 
Maguire & Steyerm. is sister to M. gattaca J. R. Grant 
from the Andes in clade III, and all Mesoamerican 
species are placed in clade II. Within the southern 


clade, M. arborescens Gilg is sister to all other species, 
which are divided into two clades: clades V and VI. 
Macrocarpaea sect. Magnolifoliae is paraphyletic in 
the sense of Grant (2005; including the viscosa clade 
of clade VI), and clade V. Clade VI corresponds to 
Macrocarpaea sect. Choriophylla. 


BIOGEOGRAPHIC ANALYSIS 


Nodes on the selected most parsimonious tree were 
classified into either allopatric (nonoverlapping), 
partially sympatric (slightly overlapping), or sympatric 
(overlapping) distributions based on the extant, 
detailed spatial distribution of each species. Past 
distributions for common ancestors to extant species 
most likely are at least somewhat different from the 
current species distribution, but because such histor- 
ical information is unavailable, no better estimate is 
possible from the current species distributions. In 
addition, the DIVA analysis provided hypotheses of 
dispersal and vicariance events that were taken into 
account when nodes were classified. The classification 
of nodes into sympatric versus allopatric species 
distributions resulted in six sympatric, 25 partially 
sympatric (— partially allopatric), and 24. allopatric 
nodes (Fig. 3). When this is done for the 19 species 
pairs in the phylogenetic tree, there are four sympatric, 
three partially sympatric, and 12 allopatric species 
pairs, given this phylogeny and species sampling. 
Allopatry is more common than true sympatry with a 
ratio of 3:1 to 4:1, but this does not account for the 
partially sympatric/allopatric distributions. 

The DIVA analysis resulted in an exact solution 
during the optimization, and an optimal reconstruction 
for our data required a total of 29 dispersal events. 
Within Macrocarpaea, there are 17 dispersal events 
and six vicariance events. Ancestral area reconstruc- 
tions and invoked dispersal events are shown in 
Figure 4. 
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species and five outgroups from Chorisepalum, Symbolan- 
thus, and Tachia.* 


Taxa/character 123456789 
Chorisepalum psychotrioides 111111000 
Symbolanthus jasonii 00-011011 
Tachia grandiflora 00-011100 
T. guianensis 00-011100 
T. occidentalis 00-012100 
Macrocarpaea angelliae 111001000 
M. apparata 111001000 
M. arborescens 00-001000 
M. auriculata 00-001000 
M. bangiana 111001000 
M. bubops 111001000 
M. chthonotropa 111001000 
M. cinchonifolia 110100000 
M. cochabambensis 111001000 
M. densiflora 222221000 
M. dies-viridis 111001000 
M. domingensis 00-001000 
M. elix 111001000 
M. ericii 00-001000 
M. fortisiana 222220000 
M. gattaca 00-001000 
M. gaudialis 00-001000 
M. glabra 00-001000 
M. glaziovii 00-001000 
M. gondoloides 111001000 
M. harlingii 111000000 
M. innarrabilis 111001000 
M. jactans 110100000 
M. jensii 111001000 
M. kuelap 111001000 
M. lenae 111001000 
M. luna-gentiana 111001000 
M. luteynii 111001000 
M. macrophylla 00-001000 
M. maguirei 110101000 
M. micrantha 111001000 
M. neblinae 00-001000 
M. nicotianifolia 00-001000 
M. noctiluca 111001000 
M. normae 22-220000 
M. obtusifolia 00-001000 
M. opulenta 111001000 
M. ostentans 22-220000 
M. pachyphylla 00-001000 
M. pachystyla 110100000 
M. pajonalis 111001000 
M. papillosa 00-001000 
M. pinetorum 00-001000 
M. pringleana 111001000 
M. revoluta 111000000 
M. robin-fosteri 110100000 
M. rubra 00-001000 
M. sodiroana 111001000 
M. stenophylla 00-001000 


Taxa/character 123456789 
M. subcaudata 111001000 
M. subsessilis 00-001000 
M. tahuantinsuyuana 110100000 
M. thamnoides 00-001000 
M. valerioi 00-001000 
M. viscosa 11001000 
M. weigendiorum 0100000 
M. zophoflora 222221090 


* Question marks indicate missing data, and hyphens 
indicate inapplicable characters. Characters and character 
states are: 1. Seed shape (3-dimensional) angular or 
spheroid (0), flattened (1). 2. Seed wings: absent (0), present 
(1). 3. Seed wings: 2-sided (0), perimetrical (1). 4. Seed 
shape (outline): square (0), long-linear (1). 5. Seed weight: 
< 0.1 mg/seed (0), > 0.1 mg/seed (1). 6. Pollen exine: 
verrucose (0), reticulate (1), smooth (2). 7. Flower position: 
terminal (0), axillary (1). 8. Corolla color: green, white, or 
yellow (0), red or purple (1). 9. Pollen aggregation: monads 
(0), tetrads or polyads (1). 


The selected optimization (Fig. 4) shows an ances- 
tral area for Macrocarpaea in southeastern Brazil (A), 
Cordillera Occidental in Colombia and Ecuador (F), 
and the Amotape—Huancabamba Zone in northern 
Peru and Ecuador (H). The first vicariance event lies 
between the Brazilian species (clade I) and Andean- 
Caribbean—Mesoamerican species. The major division 
within the Andean species into two subelades is a 
north-south vicariant event as well. The Caribbean 
species (clade IV; C) are supported as likely derived 
from Andean species (H). Dispersal patterns among 
Andean areas invoked from the DIVA analysis are 
outlined in Figure 1. Patterns within the northern 
clades (clades II and III; E, F, and G) are complicated 
and also involve one dispersal to and one from 
Mesoamerica (D), and one to the Guayana Shield (B). 
Vicariance and dispersal between the northern 
Ecuadorian and Colombian cordilleras are detected 
but represent only two dispersals: one from Cordillera 
Occidental (F) to Cordillera Central (E) and one from 
Cordillera Oriental (G) to Cordillera Central (E). 

The Amotape—Huancabamba Zone (H) is support- 
ed as the ancestral area for the southern clade. Clade 
V dispersed early southward from H to central Peru 
(I) and later to southern Peru and Bolivia (J) twice, 
with one back-dispersal northward to H. A similar 
pattern is found in clade VI, which is largely 
restricted to the Amotape-Huancabamba Zone (H) 
and shows four dispersals southward (two both to I 
and J) and one back-dispersal from J to I. Northward 
dispersals to northern Ecuador and Colombia (E and 
F) from the southern clade are only found in two 
cases. 
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Figure 1. Map of northwestern South America and southern Mesoamerica with the areas used for biogeographic analysis 
identified by letters and solid lines indicating Macrocarpaea species distributions. Arrows with numbers indicate the number 
of dispersal events as hypothesized from the DIVA analysis (see Fig. 4). Biogeographic areas are: A, southeastern Brazil; B, 
Pantepui of the Guayana Shield; C, Greater Antilles of the Caribbean; D, Mesoamerica; E, Cordillera Central in Colombia and 
Cordillera Oriental in Ecuador; F, Chocó and Cordillera Occidental in Colombia and Ecuador; G, Cordillera Oriental and 
Mérida in Colombia and Venezuela; H, Amotape—Huancabamba region in Ecuador and Peru; I, Cordillera Central in central 
Peru; J, Bolivia and Cordillera Central in southern Peru: K, Amazon Basin in Brazil, Colombia, Peru, and Venezuela. Areas A, 
C, and K are not indicated on the map. 
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Table 3. Distribution data matrix used for the DIVA 
analysis, with 0 indicating absent and 1 present for each 
specific area. See further information in Materials and 
Methods and Figure 4 for area definitions. Taxon names 
marked with * represent a larger clade: see text for 


further explanation. 


Taxon/area ABCDEFGHIJK 

Chorisepalum 01000000000 
psychotrioides 

Symbolanthus jasonii 00000001000 
Tachia grandiflora 01000000001 
T. guianensis 01000000000 
T. occidentalis 00000001111 
Macrocarpaea arborescens 00000001000 
M. auriculata* 00010000000 
M. bangiana 00000000010 
M. chthonotropa* 00000001000 
M. cinchonifolia 00000000010 
M. cochabambensis 00000000010 
M. ericii 00000001000 
M. fortisiana 00000000010 
M. gattaca 00000100000 
M. gaudialis* 00001000000 
M. glabra 00000010000 
M. gondoloides* 00000100000 
M. jactans 00000001000 
M. jensii* 00000001000 
M. luna-gentiana 00000001000 
M. macrophylla 00001100000 
M. maguirei 00000000010 
M. neblinae 01000000000 
M. nicotianifolia 00000010000 
M. noctiluca* 00000001000 
M. normae 00000000010 
M. obtusifolia* 10000000000 
M. ostentans 00000000100 
M. pachystyla 00000000100 
M. pajonalis 00000000100 
M. papillosa 00000010000 
M. pinetorum* 00100000000 
M. pringleana 00001001000 
M. revoluta 00000000100 
M. robin-fosteri 00000000100 
M. sodiroana 00000100000 
M. subsessilis* 00000001000 
M. tahuantinsuyuana 00000000100 
M. valerioi 00010000000 
M. viscosa 00000000100 
M. weigendiorum 00000000100 
M. zophoflora 00000001000 


SEEVA ANALYSIS 


Nine environmental variables were analyzed for 56 
nodes within Macrocarpaea (nodes 6 to 61 vs. the 
outgroup nodes | to 5; Fig. 3; see Table 4 for impact 
index and P values). Allopatric nodes had an average 


total impact index (i) of 0.39 (all variables and all 
nodes) whereas partly sympatric nodes had an 
average impact index of 0.26, and sympatrie nodes 
0.41. The differences in impact numbers indicate 
some differences between the groups, but there are 
only slightly larger environmental differences between 
sister clades of sympatric nodes than between sister 
clades of allopatric nodes. When terminal species 
pairs are analyzed, the results are more pronounced. 
The average impact index for allopatric species pairs 
was 0.49 (all variables and all nodes), for partly 
sympatric 0.40, and for sympatric 0.41. This finding 
indicates that allopatric species pairs may have 
ecological niches more divergent from each other 
than both sympatric and partly allopatric species. 
Allopatric terminal species pairs are the more 
ecologically different than when allopatric nodes are 
compared. 

Four Andean terminal sympatric species pairs were 
analyzed: Macrocarpaea densiflora (Benth.) Ewan 
versus M. pachyphylla Gilg in Cordillera Central of 
Colombia (Fig. 3, node 21), M. apparata J. R. Grant & 
Struwe versus M. elix J. R. Grant in southern Ecuador 
(node 37), M. bubops J. R. Grant & Struwe versus M. 
harlingii J. S. Pringle in southern Ecuador (node 38), 
and M. dies-viridis J. R. Grant versus M. lenae J. R. 
Grant in southern Ecuador (node 45). For node 21, 
only temperature seasonality showed significant 
differences between the two sister species (i — 0.54; 
P < 0.0001). Nodes 37 and 38 represent two 
sympatric species pairs that, in turn, are partially 
sympatric sister groups to each other, but the two 
species pairs show significant differences in what 
variables are different between species. Node 37 
showed significant differences in all variables except 
precipitation seasonality, with the largest differences 
in soil type (i — 1.0; P € 0.001). In contrast, node 38 
only shows significant differences in three variables: 
elevation (i — 0.34; P — 0.0037), mean annual 
temperature (i — 0.29; P — 0.0028), and minimum 


temperature of the coldest month (i — 0.25; P — 
0.0047). Node 45 has five significant variables: 
bedrock geological age (1 — 0.50; P < 0.001), 


temperature seasonality (i — 0.51; P € 0.001), annual 
precipitation (i — 0.58; P — 0.0083), precipitation of 
the driest month (i — 0.51; P € 0.0001), and 
precipitation seasonality (i — 0.51; P < 0.001). A 
similar lack of a general pattern is seen when 
allopatric or partially sympatric species pairs are 
analyzed (Table 4). 

The environmental differences between the more 
strongly supported nodes are presented here (cf. 
Fig. 2). Clades II and III are in the northern Andes, 
tepuis, and Mesoamerica, and joined by node 10 


(Fig. 3). SEEVA analysis of the two clades (Table 4) 
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Figure 2. One of 160 most parsimonious trees from the phylogenetic analysis based on combined molecular (ITS and 5S- 
NTS) and morphological data. Outgroups are Chorisepalum, Symbolanthus, and Tachia. This tree was used in the SEEVA and 
DIVA analyses. Dotted branches collapse in the strict consensus tree, and numbers below branches indicate bootstrap support 
above 50%. Current infrageneric classification for Macrocarpaea is indicated on the right. Clades marked I-VI correspond to 
nodes discussed in the text. The northern clade includes clades II + III, and the southern clade includes M. arborescens plus 
clades V + VI. 


460 


Annals of the 
Missouri Botanical Garden 


reveals that all variables show significant differences. 
Impact indices above 0.40 are found for four 
variables: elevation, soil type, annual mean temper- 
ature, and minimum temperature of coldest monih. 
Node 27 (Fig. 3) divides clades V and VI from each 
other and represents a division between a northern 
(primarily area H, clade V) and southern clade (areas 
I 4 J, clade VI). Seven of nine variables show 
differences in ecological trends between these two 
clades (Table 4), and the largest impact indices (— 
differences) are found for temperature seasonality (i — 
0.37) and precipitation seasonality (i — 0.32). Node 
42 separates a clade of three southern Ecuadorian 
species from a clade that dispersed from this area (H) 
into northern Ecuador and southern Colombia (E, F, 
and H; Macrocarpaea pringleana J. R. Grant and M. 
sodiroana Gilg). This major difference in distribution 
is associated with significant environmental differ- 
ences in six of the nine variables (67%): soil type, 
annual mean temperature, temperature seasonality, 
annual precipitation, precipitation of the driest month, 
and precipitation seasonality. The largest differences 


are found in annual precipitation (i — 0.35) and 
temperature seasonality (¢ — 0.37). 
DISCUSSION 


CONTINENTAL SCALE PATTERNS 


Similar biogeographic patterns can often be 
detected in taxa of a similar age and geographic 
origin. For South America, there are many different 
groups of this kind, from the Gondwanic relicts with 
connections to Africa, Antarctica, and Asia more than 
100 million years ago (Ma) to recent arrivals from the 
northern temperate zone along the Rocky Mountains 
of North America and the Andes when the Isthmus of 
Panama closed only a few million years ago. Other 
groups have evolved in situ on the South American 
continent and spread to North America, the Caribbean 
islands, and even farther away to New Zealand, 
Australia, and Africa. South America was, in effect, an 
isolated continent for more than 90 million years (from 
ca. 95 Ma, mid-Cretaceous Cenomanian, to ca. 4 Ma, 
Miocene; Burnham € Graham, 1999). During this 
lime, many organismal lineages were dispersing to the 
South American continent, evolved in situ, and 
dispersed within it. 

We do not know the exact age of Macrocarpaea or 
its tribe Helieae due to the lack of fossils and secure 
phylogenetic dating. Given the position of the tribe in 
the family phylogeny and an estimated age of 42— 
50 Ma for the whole Gentianaceae (Yuan et al., 2003), 
we can assume that the Helieae, and maybe the genus 
Macrocarpaea as well, has been in South America for 


at least 30 million years. Macrocarpaea is part of one 
of the more basal divergences in the tribe (the 
Macrocarpaea subclade), with only the southeast 
Brazilian genera Prepusa Mart. and Senaea Taub. 
positioned below it in the phylogeny (Struwe et al., 
2009). This minimum age estimate (30 Ma) needs to 
be tested both in a larger molecular dating analysis of 
the whole Gentianaceae family and in a detailed 
analysis of Helieae, when fossils or other additional 
evidence become available. 

The large-scale biogeographic patterns found in 
Macrocarpaea fit a scenario that is partly consistent 
with the geological history of South America. The 
DIVA analysis maps the disjunct areas (A, F, and H; 
Fig. 4) at the base of Macrocarpaea, supporting a 
disjunct ancestral area for the genus including the 
mountains of southeastern Brazil and northern Peru, 
Ecuador, and Colombia. It is too early to tell whether 
this disjunct area distribution is due to an ancient 
dispersal, or extinction in the in-between areas. The 
sister genus (Tachia) occurs in the connecting 
Amazon Basin area and is absent from southeastern 
Brazil, and therefore provides some support for a 
broad distribution of a common ancestor of Tachia and 
Macrocarpaea. An ancient dispersal event from an 
ancestor restricted to southeastern Brazil to the Andes 
is not strongly supported, because our DIVA result 
strongly supports the inclusion of the northern Andean 
areas in the ancestral area for Macrocarpaea. The 
supported scenario does not include the Bolivian 
Andes as a dispersal corridor northward from Brazil. 

The first divergence in Macrocarpaea is between 
the three species in the coastal Atlantic forest of 
southeastern Brazil, an area known for its high 
biodiversity (Myers et al., 2000), and all other species 
in the genus (Andes, tepuis, Mesoamerica, and the 
Greater Antilles). Southeastern Brazil is part of the 
ancient Gondwanic crust in South America, together 
with the Guayana Shield, and these areas represent 
relatively long-term geological stability, with the lack 
1993). 


Southeastern Brazil has repeatedly been shown to 


of historical sea incursions (Clapperton, 
include relictual, ancestrally placed lineages of 
gentians, not only in Macrocarpaea and Helieae, but 
also in Saccifolieae and Chironieae subtribe Coutou- 
beinae (Struwe et al., 2002). This southeastern Brazil— 
Andean pattern, with the Brazilian area being more 
ancestral, has also been found in other groups that 
show a large divergence in Neotropical forests (e.g., 
Fuchsia L. [Berry et al., 2004] and Gesneriaceae tribe 
Sinningieae [Perret et al., 2003, 2007)). 

Several hypotheses have been proposed for the 
derivation of the tepui flora on the Guayana Shield, 
representing long-distance dispersal from the Andes, 
derivation from lowland white-sand areas (Kubitzki, 
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1989). and in situ relictual lineages (Maguire, 1970). 
Support for several of these hypotheses is found within 
the Helieae and related tribes, with some genera (e.g.. 
Potalia Aubl.) supporting the lowland to highland 
theory (Struwe et al., 2002; Frasier et al., 2008). 

The derivation of Mesoamerican species (Macro- 
carpaea auriculata Weaver & J. R. Grant, M. 
subcaudata Ewan, and M. valerioi Standl.) from two 
lineages in the Colombian Andes is a common pattern 
in many plant groups and is likely due to dispersal 
northward along the Isthmus of Panama after its 
closing 3.1 Ma (Burnham € Graham, 1999). During 
the American interchange, many organismal groups 
moved either north or south, and other gentians that 
show a southern derivation include Potalia (Frasier et 
al., 2008) and Tachia (Struwe, unpublished data). 
According to preliminary data from DIVA-GIS 
analysis of available ecological niches (Struwe, pers. 
comm.), there are suitable habitats for Macrocarpaea 
north of Costa Rica, but these have not evidently been 
inhabited yet. 

Three species of Macrocarpaea occur on the islands 
of the Greater Antilles: M. domingensis Urb. & Ekman 
(Hispaniola), M. pinetorum Alain (Cuba), and M. 
thamnoides (Griseb.) Gilg (Jamaica), and are included 
in this study. Their biogeographic relationships are 
still uncertain because their position as being derived 
from an Andean ancestor has poor branch support in 
the phylogenetic tree, but they are strongly supported 
as being closer to Andean lineages than to the 
Brazilian lineage or Central American species. Long- 
distance dispersal from the northern Andes to the 
Caribbean is the most likely scenario. Based on the 
location of the Caribbean clade in the phylogenetic 
result, this did not happen relatively recently. The 
ancestor of the Macrocarpaea species that occurs in 
the Greater Antilles did not have wind-dispersed 
seeds, so it was most likely dispersed inadvertently by 
birds. Similar long-distance dispersal patterns be- 
tween the Andes and Hispaniola have been found 
earlier in Fuchsia (Berry, 1982, 2004). The Caribbean 
area is a complex region of several different origins, 
and its general biogeography and geological history 
are under debate and still not fully understood 
& MacPhee, 1999). Cuba and 


Hispaniola are older than Jamaica, and they were 


(Iturralde-Vinent 


connected during the early Oligocene (Iturralde- 
Vinent & MacPhee, 1999), but this was most likely 
long before the arrival of Macrocarpaea on these 
islands. 


ANDEAN PATTERNS 


The Andean uplift started in the Miocene (ca. 
20 Ma) and continued until the Holocene, generally 


moving from the southern part of the mountain range 
toward the northern part, with the main northern 
upheaval ca. 2-8 Ma (Haffer, 1987). Other scientists 
support an earlier start at 40 Ma, with major northern 
uplift at 18 Ma, as suggested by Ghosh et al. (2006) 
and Gregory-Wodzicki (2000). The Andean forests of 
the northern Andes where Macrocarpaea occurs have 
been hypothesized to be of Miocene to lower Pliocene 
in age (van der Hammen, 1979), but their range in 
elevation shifted downward during the Pleistocene 
glaciations (Haffer, 1987). The Andean clade of 
Macrocarpaea is clearly separated into two subclades, 
a northern and southern one, divided at the Amotape— 
Huancabamba Zone. The geographic division of these 
two groups follows the Amotape Cross, a geological 
shear zone, which divides the Andes (and Ecuador) 
into two parts: the younger northern Andes (formed in 
Late Pliocene to Pleistocene) and older central Andes 
(Miocene to Pliocene; Young & Reynel, 1997). 

Based on our DIVA analysis, the ancestral 
distribution within the Andes includes the Amo- 
tape-Huancabamba Zone (H), which is among the 
most species-rich areas in the Andes. Also included in 
the ancestral Andean area are the Cordillera Occi- 
dental of Colombia and Ecuador (F) and the Cordillera 
Oriental of Colombia and Cordillera de Mérida in 
Venezuela (G). The deep split between the northern 
and southern clades of Macrocarpaea represents an 
ancient divergence that strongly correlates with 
current distribution patterns. The fact that only two 
dispersal events have crossed over from the south to 
the north, and none from the north to the south, is 
remarkable. This zone in central and northern 
Ecuador is also known for its absence of Tachia and 
Symbolanthus, two other Helieae genera that have 
suitable habitats in the boundary area. Possible 
explanations for this deep divide might be the early 
division of the Andes in this region between a 
northern and southern part, since they have been 
divided by sea incursions (early Miocene), mountain 
uplift and creation of valleys, and/or volcanism 
(Clapperton, 1993; Burnham & Graham, 1999), or 
perhaps a westward flow of the Amazon into the 
Pacific Ocean (Mapes et al., 2006). 

Quaternary volcanism is absent from the Amotape— 
Huancabamba Zone (area H) and the central Peruvian 
area (H and J), but occurs in areas both north and 
south of these regions (Clapperton, 1993). Volcanism 
could affect extinction rates severely. and the lack of 
voleanism might partly explain the relatively higher 
species numbers in the Amotape—Huancabamba 
Zone. In addition, this area is characterized by humid 
isolated forest islands occurring in a highly dissected 
landscape (Jérgensen, pers. comm.). The species of 
the southern clade also have wind-dispersed, winged 
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Figure 3. The selected most parsimonious tree used in DIVA and SEEVA analyses, with numbered nodes and geographic 
distribution listed for each species. Biogeographic area coding for DIVA is indicated with letters A-K after species names (see 


text for definition of areas). Symbols (O, 


. %) on nodes indicate sympatric, partially sympatric, or allopatric clades, based 


on exact species distribution within Macrocarpaea (not DIVA area classification), which is further analyzed in the 


SEEVA analysis. 
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Figure 4. Results of the DIVA analysis when mapped onto the selected most parsimonious phylogenetic tree. Areas 
optimized onto each branch are marked with letters (see legend and Fig. 3 for coding): when several optimizations were equally 
optimal, one was selected based on spatial and geological information (see Appendix 1 for all optimal DIVA reconstructions). 
Circles (B) on branches indicate dispersal events leading to speciation or enlargement of geographic distribution for a species, 
and X by a node indicates a vicariance event leading to speciation. The schematic diagram to the left shows an overview of area 
relationships, with bold letters (D, G) indicating expansion to new areas, solid lines indicating vicariance events, and arrows 
indicating dispersal events. Vicariances and dispersals are only marked for Macrocarpaea. 
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seeds that presumably could lead to increased 
colonization of new areas and thereby increased 
allopatric speciation through island hopping. Howev- 
er, more instances of allopatric speciation are found in 
the northern clade than in the southern clade. 
Weigend (2002, 2004) showed that the Huancabamba 
depression is not a dispersal barrier to Andean 
species of middle elevations, and this conclusion is 
upheld by our data. This finding is in contrast to data 
from Eleutherodactylus Duméril & Bibron frogs 
(Young & Reynel, 1997; Duellman & Pramuk, 
1999); however, it should be noted that this vertebrate 
group represents an immense Neotropical radiation 
with a very different speciation history. 


SPECIATION PATTERNS 


Within Andean Macrocarpaea, allopatric speciation 
is much more common than sympatric speciation, and 
this has also been shown for Gesneriaceae in 
southeastern Brazil (Perret et al., 2007). In some 
cases, allopatric speciation is linked to long-distance 
dispersal (e.g., M. neblinae on Sierra de Neblina in 
southern Venezuela, and M. valerioi in Costa Rica), 
but often allopatric sister taxa occur in more adjacent 
parts of the same mountain range (e.g., M. cinchoni- 
folia (Gilg) Weaver in western Bolivia and M. normae 
J. R. Grant in southern Peru). Such north-south 
patterns are common and indicate either dispersals 
along mountain ridges or historical splitting of larger 
ancestral populations into separate species after 
isolation. Both scenarios support the theory that 
mountain ranges such as the Andes represent virtual 
islands (Young et al., 2002). In 


Macrocarpaea, most species occur on the slopes and 


the case of 


not on the summits, but slopes may serve as isolated 
units as well. Isolated páramos on Andean summits 
have long been considered analogs to island chains 
such as Hawaii and the Galápagos (Young et al., 
2002}, but our data also support this for lower 
elevations. Many of the paramo plants are relatively 
recent immigrants from plant groups from the northern 
temperate zones, whereas the forested slopes include 
primarily Neotropical floristic elements that reach 
high species diversity in the Andes. Within Helieae, 
the highest species diversity is found in the Andes, 
but is largely represented by only two genera, 
Macrocarpaea and Symbolanthus (Molina & Struwe, 
2008). The Amazon lowlands and the Brazilian and 
Guayana Shields have much fewer species, but they 
represent the ancestral evolutionary lineages (Struwe 
et al., 2002, 2009). 

The seeds of Macrocarpaea can be divided into two 
major types: very small and winged or larger, angular, 
and heavier. Winged seeds, which promote dispersal 


over larger distances, are found in clades V and VI, 
which show repeated dispersal patterns north-south 
within the southern clade. Such seeds are also found 
in the northern group at node 12 (Fig. 3), M. 
gondoloides J. R. Grant and M. luteynii J. R. Grant 
& Struwe, as an independently derived character trait. 
The northern group with heavier seeds shows fewer 
occasions of dispersal, but is also more poorly 
sampled. Because it was not possible to include all 
species in our study due to unavailability of material, 
species from the southern area are overrepresented. 
Pollinator information for Macrocarpaea species is 
relatively scarce, but most species appear to be 
generalists and are pollinated by a multitude of 
animal groups and species (Grant, pers. obs.) 
Pollinator segregation therefore does not support 
sympatric speciation in Macrocarpaea, and it is more 
likely that ecological niche divergence or population 
isolation and subsequent fixation of different traits in 
smaller populations have led to different species. 


ECOLOGICAL NICHE PATTERNS 


Results from the SEEVA analysis of sister species 
show that all species pairs have individual divergence 
patterns and environmental differences. Divergence in 
ecological niches is common both in allopatric and 
sympatric species in Macrocarpaea, and generally 
these divergences are not significantly different in 
size between the two speciation types, but these 
variables differ between species pairs. For example, 
divergences were found based on different altitudinal 
zones (M. apparata and M. elix; node 37, Fig. 3) and 
on different types of bedrock and in zones with 
different climate seasonality (M. dies-viridis and M. 
lenae; node 45, Fig. 3). Perret et al. (2007) also found 
a lack of increased divergence in sympatric species 
when compared with allopatric species. 

The species differences found with SEEVA show 
only patterns, not processes, and variables should not 
be seen as the probable cause for speciation, unless 
further studies can show adaptation to specific 
environments or changes linked to paleoclimatological 
or geological events. We know that species have 
moved around on a geographic scale, especially 
during the Pleistocene, but their ecological niches 
might have been more stable due to niche conserva- 
tism. Using a different approach, Peterson et al. 
(1999) showed that speciation through geographic 
separation often appears before ecological niche 
separation in vertebrates and butterflies in Mexico, 
and our results indicate the same. Extracting 
environmental data from current locations therefore 
most likely represents historical ecological niches, if 
not the historical location of the population. 
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In future studies, it would be interesting to compare 
species pairs for particular areas from many different 
genera to determine whether ancestral populations 
reacted similarly to climatological and geological 
events in the past, and whether recurrent ecological 
niche divergence appeared in unrelated lineages. This 
would be possible, for example, in the Amotape— 
Huancabamba Zone, where we now have at least four 
different data sets from angiosperms (including the 
studies by Weigend, 2002, 2004). 

Our study shows that explanations of speciation 
events need to be sought individually for each species 
that 
applicable across a larger species group distributed 


pair, and generalities are not necessarily 


over a large area. If we want to understand speciation 
patterns, the ecological niches of species, and threats 
and means to the conservation of these species, much 
more data need to be collected and analyzed (Young 
et al., 2002). One major difficulty is the lack of up-to- 
date revisionary, phylogenetic, and georeferenced 
data for most plant genera. Species-level phylogenies 
only represent a small percentage of Neotropical 
biodiversity, and there is a dire need for more 
taxonomic work that can be integrated with biogeog- 
raphy, ecology, and conservation. 
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APPENDIX ] 


Optimal distributions from the DIVA analysis listed for 
each node (ancestor of terminals [AOT]; node numbers in 
Fig. 3: letters refer to legend for Fig. 4). Note that 
optimizations in Figure 4 sometimes reflect a dispersal event 
at a terminal branch that might be indicated as a joint 
distribution in the list below, because changes in area 
distributions on terminal branches are not mapped with 
DIVA. Optimizations outside of Macrocarpaea are not listed 
here. and due to the collapse of groups with homogenous 
areas (see Materials and Methods). some node numbers are 
not listed here. 

Node 6. AOT M. obtusifolia-M. normae: AB AC ABC AD 
ABD ACD ABE AF ABF ACF ADF AG ABG ACG ADG 
AFG AH ABH ACH ADH AFH AGH 

Node 9, AOT M. glabra-M. normae: CD CF CDF CG CDG 
CFG DH CDH FH CFH DFH GH CGH DGH FGH 

Node 10, AOT M. glabra-M. macrophylla: F G DG FG DFG 
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ode 11, AOT M. gondoloides-M. macrophylla: F DF FG 
DFG EFG 

ode 13, AOT M. valerioi-M. macrophylla: D DF DEF DG 
EG DEG FG DFG EFG 

ode 14, AOT M. auriculata-M. macrophylla: DE DF DEF 
ode 16, AOT M. valerioi-M. papillosa: DG 

ode 17, AOT M. glabra—M. nicotianifolia: G FG BFG EFG 
ode 18, AOT M. neblinae-M. nicotianifolia: EF BEF BG 
BEG FG BFG EFG 

ode 19, AOT M. gaudialis-M. nicotianifolia: EG 

ode 22, AOT M. neblinae-M. gattaca: BF 

ode 23, AOT M. pinetorum—M. normae: CH 

ode 26, AOT M. arborescens—M. normae: H 

ode 27, AOT M. viscosa—M. normae: HI HJ HIJ 

ode 28, AOT M. viscosa—M. pajonalis: H 

ode 29, AOT M. bangiana—M. pajonalis: HJ 

Node 30, AOT M. subsessilis-M. pajonalis: H 

Node 32, AOT M. jensii-M. pajonalis: H 


Node 33, 
Node 34, 
Node 39, 
Node 43, 
Node 44, 
Node 49, 
Node 50, 
Node 51, 
Node 52, 
Node 53, 
Node 54, 
Node 55, 
Node 56, 
Node 57, 
Node 58, 
Node 59, 
Node 60, 
Node 61, 


AOT M. 
AOT M. 
AOT M. 
AOT M. 


chthonotropa—M. pajonalis: H 
noctiluca—M. pajonalis: H 
ericii-M. pajonalis: HI 
pringleana—M. sodiroana: FH EFH 


AOT M. jensii-M. sodiroana: H 


AOT M. 
AOT M. 
AOT M. 


viscosa-M. pachystyla: HI 
luna-gentiana-M. pachystyla: H 
zophoflora—M. pachystyla: HJ 


AOT M. fortisiana—M. pachystyla: J 


AOT M. 
AOT M. 
AOT M. 
AOT M. 
AOT M. 
AOT M. 
AOT M. 
AOT M. 
AOT M. 


maguirei-M. pachystyla: IJ 
cochabambensis-M. normae: IJ 
revoluta-M. normae: | 
weigendiorum—M. normae: | 
ostentans-M. normae: 1 IJ 
robin-fosteri-M. normae: | 
robin-fosteri-M. jactans: HI 
tahuantinsuyuana—M. normae: IJ 
cinchonifolia-M. normae: J 


